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摘要 : 本 文采 用 分 子 动力 学 模拟 方法 研究 离子 的 存在 对 水 蒸气 核 化 过 程 的 影响 。 模 拟 系 综 初 始 温度 为 600 K， 达 到 平衡 后 撤去 
水 蒸气 热 滩 ， 将 载 气温 度 标 定 至 200 K， 与 水 莹 气 进行 热量 交换 ， 带 走 水 蒸气 核 化 释放 的 核 化 潜 热 。 研 究 结果 表明 ， 离 子 的 存 
在 能 够 促使 核 化 的 发 生 ， 但 离子 的 浓度 对 系 综 内 最 大 团 簇 的 生长 影响 较 小 。 根 据 团 簇 分 析 统 计 系 综 内 含 20 个 NaVCL 时 ， 核 化 
率 约 为 1.01x1028/(cms.s)， 与 经 典 核 化 理论 计算 值 4.89x1028/(cms.s) 在 同一 个 数量 级 。 

关键 词 : 离子 ， 蜡 质 核 化 ， 团 簇 ; 分子 动力 学 

中 图 分 类 号 : TK123 文献 标识 码 : A 


Molecular Dynamic Study of the Effect of Ions on Water Vapor Nucleation 
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Abstract: Molecular dynamics Simulation was employed to explore the mechanism of heterogeneous nucleation of 
water vapor under the effect of ions. After the System reached to the equilibrium state at initial temperature 600 K, 
the temperature of carrier gas was rescale at 200 K, which worked as the heat resource for the water vapor nucleation. 
The nucleation heat was removed by carrier gases during the nucleation process. The results show that, ions promote 
the nucleation progress on the initial stage, and the concentration of ions has less effect on the growth of the largest 
cluster in the system. The nucleation rate obtained from the cluster analysis was 1.01x1023/(cm3.s) in the system with 
20 Nat and 20 CL ions, which is in the same order of magnitude with the value 4.89x1028/(cm3.s) predicted by 
classic nucleation theory analysis. 
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子 是 几 个 中 性 分 子 团聚 在 一 个 带电 离子 周围 形成 的 
粒子 ， 小 粒子 吸附 在 较 大 的 中 性 气 溶胶 粒子 上 形成 
相 变 是 自然 界 广泛 存在 的 一 种 物理 现象 ， 在 能 ”大 离子 。 同 一 符号 的 大 小 离子 的 浓度 在 海洋 乡村 城 
源 利用 、 强 化 传 热 等 领域 起 着 重要 作用 023。 在 大 气 。 市 不 尽 相 同 ， 海 洋 上 空 小 离子 较 多 ， 城 市 上 空 大 离 
层 中 ， 当 水 蒸气 的 相对 饱和 度 达 到 过 饱和 或 大 气 中 ” 子 较 多 ， 并 且 大 小 离子 浓度 变化 趋势 相反 。 海 陆 上 
有 足够 多 的 灰尘 粒子 时 ， 就 会 发 生 水 蒸气 凝结 过 程 。 空 的 正 负 离子 浓度 相差 不 大 ， 但 正 离子 稍 多 B]。 
而 形成 云雾 。 随 着 经 济 的 发 展 、 科 技 水 平 的 提高 ， Winkler 等 由 实验 研究 有 机 蒸气 在 大 小 为 1-24 
气 溶胶 颗粒 物 污染 越 来 越 引起 人 们 的 注意 ， 大 气 气 。 nm 的 小 团 复 和 纳米 颗粒 上 的 异 质 核 化 过 程 并 与 经 
溶胶 严重 的 影响 了 空气 质量 、 危 害 了 人 类 健康 。 大 ” 典 核 化 理论 及 Kelvin 关系 对 比 ， 结 果 发 现 ， 异 质 核 
气 中 悬浮 的 离子 义 分 为 小 离子 与 大 离子 两 种 ， 小 离 。 越 小 越 能 激发 核 化 。Inci 等 5 通过 引入 与 核 化 蒸气 分 
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子 大 小 相同 的 粒子 作为 异 质 核 来 研究 异 质 核 化 ， 并 
改变 异 质 核 与 蒸气 分 子 的 相互 作用 力 来 讨论 其 对 核 
化 的 影响 。 研 究 结 果 表明 ， 异 质 核 化 与 蒸气 分 子 相 
互 作用 力 较 小 的 时 候 ， 异 质 核对 小 团 簇 胚 核 的 形成 
功 影响 很 小 ， 随 着 相互 作用 力 的 增加 ， 团 徐 形 成 自 
由 能 降低 。 张 超 等 1 采用 团 簇 分 析 法 分 析 大 于 不 同 阐 
值 的 团 复数 随时 间 的 变化 情况 ， 通 过 最 小 二 乘法 分 
别 拟 合 获取 3 种 模拟 工 况 下 的 核 化 速率 ， 对 比 不 同 
带电 量 离子 对 核 化 过 程 的 影响 。 结 果 表 明正 二 价 离 
子 比 正 一 价 离子 周围 产生 的 电场 力 大 ， 对 核 化 过 程 
影响 大 。 因 此 离子 的 存在 对 水 蒸气 核 化 有 着 重要 的 
影响 ， 而 正 负 离子 及 离子 浓度 对 水 蒸气 核 化 的 微观 
机 理 还 有 待 于 进一步 的 研究 ， 本 文 建立 分 子 动力 学 
模拟 模型 ， 研 究 离 子 影 响 水 蒸气 核 化 的 微观 机 理 及 
核 化 的 微观 物理 过 程 。 


2 模拟 模型 与 计算 细节 


采用 Nar+、CL 离 子 代 表 正 负 粒 子 研 究 其 对 水 蒸气 
核 化 的 影响 。 图 1 为 带电 离子 存在 下 的 水 蒸气 核 化 
模拟 系统 图 ，5000 个 SPC/E 模型 的 水 蒸气 分 子 均匀 
分 布 在 模拟 盒 29x29x29 nm3 内 ，5000 个 氮 蒸 气 分 子 
作为 载 气 (未 画 出 )， 系 统 平衡 后 控制 载 气 温度 ， 载 
气 可 带 走 水 蒸气 核 化 释放 的 显 热 和 潜 热 。 为 使 图 示 
清晰 ， 水 分 子 用 点 表示 ，Na+、CIL 离 子 用 放大 球体 表 


修 。 


图 1 存在 带电 离子 的 水 蒸气 核 化 模型 


Fig. 1 Schematic plot of water nucleation System with ions 


氯 分 子 间 相互 作用 采用 经 典 的 L-J (12-6) 势 函 
数 。 水 分 子 选 用 SPC/E 模型 ， 作 用 势 函 数 为 ， 
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由 于 氧 原 子 质量 较 小 ， 方 程 (1) 中 第 一 项 L-J 势 
函数 只 用 来 计算 氧 原 子 间 的 相互 作用 力 。 Ooo 和 coo 
分 别 为 氧 原子 间作 用 的 尺度 参数 和 势 阱 深度 。 后 一 
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项 计算 各 个 带电 荷 原子 间 的 静电 力 ,9 与 人 分 别 为 原 
子 i 与 j 的 电荷 , rj 为 原子 i 与 j 间 的 距离 。go 为 真空 
介 电 常数 。 水 分 子 与 握 原 子 间 的 相互 作用 力 由 修正 
的 L-J 势 函 数 计算 ， 水 分 子 与 Nat+、CI 离子 的 作用 
势 函 数 取 自 文献 中 ,根据 Lorentz-Berthelot 混合 法 则 
计算 其 与 水 分 子 上 的 氧 原子 作用 的 能 量 参数 与 尺度 
参数 。 表 1 给 出 势 函 数 中 各 原子 作用 参数 及 电 葆 。 
时 间 步 长 3 fs, 截断 半径 1.5 nm, 采用 PPPM 方法 对 
长 程 静电 力 进行 修正 ,采用 周期 性 边界 条 件 , Velocity 
Verlet 算法 求解 分 子 动力 学 方程 。 首 先 ， 采 用 
Nose-Hoover 热 浴 , 模拟 300 ps 使 得 系统 在 600K 温 
度 下 达到 平衡 ， 撤 去 热 浴 将 载 气温 度 标定 至 200 K。 
再 模拟 2000 ps 统计 参数 分 析 核 化 过 程 。 


表 1. 势 函数 中 各 原子 作用 参数 及 电荷 


Table 1 Parameters and charges of atoms for potential 


function 
Atom /人 A 2 /kJ mol! gile 
O 3.166 0.6502 -0.8476 
H 0 0 0.4238 
Ar 3.405 0.996 0.0 
Na 2.583 0.4184 1.0 
Cl 4.400 0.4184 -1.0 


3 模拟 结果 及 分 析 
3. 1 核 化 过 程 热力 学 参数 

为 讨论 离子 浓度 对 水 蒸气 核 化 的 影响 ， 系 综 内 
分 别 引 入 20、50、100、200 个 NarCL 离 子 对 。 图 2 
给 出 了 不 同 离子 浓度 系 综 内 水 蒸气 的 动能 与 势能 随 
时 间 的 变化 关系 。 系 综 在 模拟 达到 平衡 后 撤去 热 浴 ， 
将 载 气温 度 标定 至 200 K， 在 载 气 热 浴 下 ， 水 蒸气 温 
度 降 低 至 核 化 平衡 温度 。 由 图 2(a) 温 度 的 分 布 图 可 以 
看 出 ， 在 核 化 初期 ， 离 子 浓度 较 低 的 系统 ， 核 化 慢 ， 
释放 潜 热 少 ， 温 度 下 降 略 快 ; 离子 浓度 较 高 的 系统 ， 
系 综 内 水 蒸气 核 化 较 剧烈 ， 核 化 释放 的 潜 热 未 能 被 
及 时 带 走 ， 温 度 降 低 略 慢 。 随 着 核 化 进行 ， 温 度 均 
趋 于 平衡 值 ， 系 综 内 离子 浓度 对 平衡 温度 的 影响 较 
小 ， 其 对 应 的 核 化 平衡 温度 为 372~375 K。 核 化 过 程 
是 势能 降低 的 过 程 ， 由 于 离子 的 相互 吸引 使 得 分 子 
间距 离 减 小 ， 使 得 离子 浓度 较 大 的 系 综 对 应 的 势能 
较 低 。 
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图 2 离子 存在 时 水 蒸气 核 化 过 程 热力 学 参数 随时 间 变 化 关系 


Fig. 2 Evolution of thermal properties in nucleation of water 


vapor with ions (a) temperature (b) potential energy 
3.2 核 化 过 程 分 析 

Na+、CI 与 水 分 子 间 的 短程 范 德 华 力 相对 较 弱 ， 
其 主要 的 相互 作用 力 来 之 长 程 库伦 力 。 图 3 为 系 综 
内 包含 100NaVCT 正 负离子 时 水 燕 气 核 化 的 动态 过 
程 。 与 之 对 应 图 4 为 系统 内 不 同时 刻 点 团 簇 大 小 的 
分 布 。 系 统 模拟 300 ps 使 其 达到 600 K 高 温 的 平衡 
状态 ， 此 时 水 落 气 分 子 均匀 地 分 布 在 系统 内 ， 无 稳 
定 团 簇 存在 。 然 后 撤去 热 浴 ， 将 载 气 温度 标定 至 200 
K。 随 着 载 气 与 水 蒸气 进行 热量 交换 ， 带 走 显 热 ， 水 
蒸气 系统 温度 降低 。 核 化 为 亚 稳 态 过 程 ， 离 子 与 水 
分 子 的 较 强 作用 力 相当 于 一 些 扰动 ， 促 使 核 化 发 生 ， 
温度 降低 到 核 化 点 时 ， 在 离子 周围 的 水 分 子 以 离子 
为 核 化 中 心 核 化 。 由 于 原子 间 的 相互 作用 随 着 距离 
的 增加 而 减弱 ， 所 以 带电 离子 对 其 周围 水 分 子 核 化 
的 影响 也 随 距 离 变 化 。 如 图 3 中 t=500 ps 时 水 蒸气 
分 子 核 化 形成 团 复 的 情况 , 对 应 团 徐 大 小 分 布 见 图 4 
中 t=500 ps， 这 时 系统 内 大 量 小 团 复 形成 ， 不 稳定 的 
小 团 复 在 分 子 相互 作用 下 破裂 ， 稳 定 的 小 团 簇 继续 
长 大 。 当 tf1300 ps， 形 成 的 小 团 徐 开始 在 相互 作用 
下 聚 并 形成 大 团 徐 ， 最 大 团 复 约 由 500 个 水 分 子 组 
成 ， 此 时 仍 有 一 些小 的 不 稳定 小 团 复 存在， 不 断 地 
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破裂 又 形成 , 系统 内 水 燕 气 饱 和 度 降 低 。 当 t=2000 ps 
时 ， 系 统 内 的 核 化 基本 达到 一 个 动态 平衡 ， 几 个 较 
大 的 团 簇 稳 定 存 在 ， 最 大 团 簇 约 包含 700 个 水 分 子 ， 
水 分 子 单 体 均匀 分 布 在 系统 内 。 继续 模拟 到 2300 ps， 
系统 内 的 团 簇 大 小 分 布 变 化 较 小 ， 说 明 核 化 系统 达 
到 了 核 化 与 蒸发 的 动态 平衡 状态 。 


t=2000 ps do 


t=2300 ps 


图 3 系 综 内 水 蒸气 核 化 的 动态 过 程 〈100NarWCL 离 子 ) 


Fig. 3 The snapshots for dynamic process of water vapor 


nucleation with 100 Nat ions and 100Cl ions 


T T T 3 


t=500 ps 


T T T 3 


{=1300 ps 


200 400 600 800 
T T T 


{=2300 ps 


Cluster size distribution 


HD 
400 600 800 


Cluster size 


图 4 系 综 内 团 簇 大 小 随时 间 的 分 布 (100Nat*/CIl 离子) 


Fig. 4 Evolution of cluster size distribution in water nucleation 


process with 100 Nar ions and 100Cl ions 
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图 5 系统 内 最 大 团 簇 的 大 小 随时 间 的 分 布 


Fig. $ Evolution of the largest cluster size distribution in water 


nucleation process with ions 

同时 讨论 了 离子 浓度 对 水 蒸气 核 化 的 影响 ， 分 
析 不 同 离子 浓度 下 ， 系 统 内 最 大 团 徐 包含 的 分 子 数 
随 核 化 时 间 的 变化 ， 如 图 5 所 示 。 在 核 化 初期 ， 离 
子 的 存在 能 够 促使 核 化 的 发 生 ， 团 簇生 长 较 快 ， 图 
发 生 阶 路 为 小 团 簇 聚 并 到 大 团 禾 ， 最 大 团 驴 急剧 
长 大 。 随 着 核 化 进行 ， 小 团 簇 聚 并 长 大 ， 系 统 达 到 
新 的 平衡 状态 ， 最 大 团 艇 包含 约 700 水 分 子 。 不 同 
离子 浓度 的 系统 ， 核 化 达到 平衡 时 ， 最 大 团 徐 大 小 
基本 相同 ， 离 子 浓度 影响 较 小 。 
3. 3 核 化 率 计算 

核 化 率 定义 为 单位 时 间 单 位 体积 内 生产 的 活化 
核 ( 即 能 够 稳定 存在 并 长 大 的 分 子 团 艇 )。 本 文 以 包含 
20Nat/CIL 离子 的 核 化 系 综 为 例 ， 分 别 采用 经 典 核 化 
理论 预测 及 团 簇 分 布 统计 分 析 离 子 存在 时 水 蒸气 的 
核 化 率 大 小 。 核 化 系统 内 对 应 的 饱和 度 约 为 2.9。 离 
子 存 在 时 水 蒸气 核 化 达到 动态 平衡 时 对 应 的 温度 约 
为 375 K, 此 温度 下 对 应 的 饱和 气 与 饱和 液 的 密度 分 
别 为 6.27x104 gcm3 和 0.9165 gcm-3, 参考 表面 张力 
36.94 Nm 。 由 于 本 研究 中 的 带电 离子 没有 考虑 大 
小 ， 采 用 均 质 核 化 的 核 化 率 计 算 公 式 ， 


2 AG, i 
Ja = PrP’ exp(— oe) (2) 
B 


计算 得 出 的 核 化 率 的 大 小 为 4.89x1023/(cm3.s)。 
此 外 根据 核 化 率 定义 ， 可 从 团 簇 分 析 方 法 求 得 核 化 
率 ， 图 6 给 出 了 包含 分 子 数 大 于 一 定 阔 值 的 团 禾 团 
复数 随时 间 的 变化 ， 阅 值 分 别 为 10、20、30、40、 
50。 当 阔 值 为 10 时 ， 系 综 内 团 簇 数 先 迅 速 增加 后 降 
低 ， 变 化 不 稳定 。 当 阔 值 大 于 20 时 ， 团 徐 数 增加 的 
斜率 基本 相同 ， 从 而 可 以 根据 其 变化 的 斜率 统计 获 
得 核 化 率 。 图 中 的 虚线 为 团 徐 生长 阶段 的 线性 拟 合 ， 
此 时 对 应 的 系 综 内 包含 20NarCL 离 子 时 水 蒸气 核 化 
的 核 化 率 约 为 1.01x1028Kcms s)， 与 经 典 核 化 理论 预 
测 值 在 同一 数量 级 , 高 于 文献 四 计算 得 出 的 均 质 核 化 
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的 核 化 率 值 2.8x102 /cms.s)， 说 明 离 子 的 存在 对 水 
蒸气 核 化 有 促进 作用 。 采 用 同样 的 方法 根据 团 徐 分 
析 得 出 系 综 内 包含 50、100、200 NarCL 离 子 时 对 应 
的 核 化 率 依次 为 ，1.54x1023/(cm3 s)、2.27x1023/(cm3 
s)、4.52x1023/(cm3 s)， 核 化 率 随 离子 浓度 的 增加 而 非 
线性 的 增加 。 
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图 6 系 综 内 团 簇 数 随时 间 的 变化 (20Nat/CI) 


Fig. 6 Evolution of clusters larger than a threshold size in the 


nucleation process with 20 Nar and 20CL ions 
4 结论 


本 文 建立 了 Nat/Cl 离 子 存在 时 水 蒸气 核 化 的 分 
子 动力 学 模型 ， 模 拟 研究 Nat/Cl 离 子 及 离子 浓度 对 
水 燕 气 核 化 的 影响 ， 分 析 核 化 过 程 温 度 、 势 能 的 变 
化 ， 根 据 团 簇 分 析 及 经 典 核 化 理论 预测 离子 存在 时 
水 蒸气 的 核 化 率 。 结 果 发 现 ， 在 核 化 初期 ， 离 子 的 
存在 能 够 促使 核 化 发 生 ， 但 随 着 核 化 进行 至 动态 平 
衡 时 ， 离 子 浓 度 对 系 综 内 最 大 团 复 大 小 影响 较 小 。 
以 含 20Nat/CL 离 子 系统 为 例 ， 根 据 核 化 率 定义 由 团 
簇 分 析 统 计 得 出 核 化 率 约 为 1.01x1028(cms s)， 高 于 
均 质 核 化 系统 ， 且 与 经 典 核 化 理论 计算 值 
4.89x1028/(cms s) 在 同一 个 数量 级 。 
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